home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 1998 November / IRIX 6.5.2 Base Documentation November 1998.img / usr / share / catman / p_man / cat3 / complib / CHPFA.z / CHPFA
Text File  |  1998-10-30  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCHHHHPPPPFFFFAAAA((((3333FFFF))))                                                            CCCCHHHHPPPPFFFFAAAA((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CHPFA   - CHPFA factors a complex Hermitian matrix stored in packed form
  10.      by elimination with symmetric pivoting.
  11.  
  12.      To solve  A*X = B , follow CHPFA by CHPSL.  To compute  INVERSE(A)*C ,
  13.      follow CHPFA by CHPSL.  To compute  DETERMINANT(A) , follow CHPFA by
  14.      CHPDI.  To compute  INERTIA(A) , follow CHPFA by CHPDI.  To compute
  15.      INVERSE(A) , follow CHPFA by CHPDI.
  16.  
  17.  
  18. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  19.       SUBROUTINE CHPFA(AP,N,KPVT,INFO)
  20.  
  21. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  22.      On Entry
  23.  
  24.      AAAAPPPP COMPLEX (N*(N+1)/2)
  25.         the packed form of a Hermitian matrix  A .  The
  26.         columns of the upper triangle are stored sequentially
  27.         in a one-dimensional array of length  N*(N+1)/2 .
  28.         See comments below for details.
  29.  
  30.      NNNN INTEGER
  31.         the order of the matrix  A .  Output
  32.  
  33.      AAAAPPPP A block diagonal matrix and the multipliers which
  34.         were used to obtain it stored in packed form.
  35.         The factorization can be written  A = U*D*CTRANS(U)
  36.         where  U  is a product of permutation and unit
  37.         upper triangular matrices , CTRANS(U) is the
  38.         conjugate transpose of  U , and  D  is block diagonal
  39.         with 1 by 1 and 2 by 2 blocks.  KVPT    INTEGER(N)
  40.         an integer vector of pivot indices.
  41.  
  42.      IIIINNNNFFFFOOOO INTEGER
  43.         = 0  normal value.
  44.         = K  if the K-th pivot block is singular.  This is
  45.         not an error condition for this subroutine,
  46.         but it does indicate that CHPSL or CHPDI may
  47.         divide by zero if called.  Packed Storage The following program
  48.      segment will pack the upper triangle of a Hermitian matrix.
  49.         K = 0
  50.         DO 20 J = 1, N
  51.         DO 10 I = 1, J
  52.         K = K + 1
  53.         AP(K)  = A(I,J)
  54.         10    CONTINUE
  55.         20 CONTINUE LINPACK.  This version dated 08/14/78 .  James Bunch,
  56.      Univ. Calif. San Diego, Argonne Nat. Lab.  Subroutines and Functions BLAS
  57.      CAXPY,CSWAP,ICAMAX Fortran ABS,AIMAG,AMAX1,CMPLX,CONJG,REAL,SQRT
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.